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For  an  integer  program,  a  k-cut  is  a  cutting  plane  generated  by  the  Gomory 
mixed  integer  procedure  from  a  row  of  the  LP  tableau  after  multiplying  it  by 
a  positive  integer  k.  With  this  terminology,  Gomory  mixed  integer  cuts  are 
just  1-cuts.  In  this  paper,  we  compare  the  &:-cuts  (k  >  2)  with  Gomory  mixed 
integer  cuts.  In  particular,  we  prove  in  the  pure  case  that  with  exactly  50% 
probability  the  k- cuts  perform  better  variable-wise  than  the  Gomory  mixed 
integer  cuts,  and  vice  versa.  Some  computational  experiments  on  knapsack 
problems  are  reported  to  illustrate  this  property. 

(. Integer  Programming;  LP  Tableau;  Gomory  Mixed  Integer  Cut;  K-cut ) 


1.  Introduction 

In  the  late  fifties  and  early  sixties,  Gomory  introduced  two  kinds  of  general 
cutting  planes,  fractional  cuts  and  mixed  integer  cuts  [9,  10,  11],  for  solving 
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integer  programming  problems.  Later  on,  extensive  work  was  done  on  the 
properties  of  Gornory  cuts  and  the  convergence  of  Gomory  cutting  plane 
methods  [3,  5,  7,  8,  14,  22,  23].  The  potential  of  Gomory  cuts  in  computation 
was  also  explored  [20,  21].  The  advantage  of  Gomory  cuts  within  a  branch- 
and-cut  framework  was  computationally  shown  in  [2], 

Glover  [8]  proposed  to  linearly  combine  the  original  equations  from  the 
linear  programing  tableau  with  Gomory  cuts.  The  combination  allows  two 
free  parameters  to  be  fixed.  Glover  discussed  how  to  specify  the  parame¬ 
ter  values  so  that  the  resulting  cuts  “not  only  limit  the  feasible  set  more 
restrictively  than  the  method  of  integer  form,  but  satisfy  other  criteria  as 
well”. 

In  a  series  of  papers  [12,  15,  16,  18]  in  the  late  sixties  and  early  seventies, 
Gomory  and  Johnson  showed  how  to  generate  cutting  planes  from  subaddi¬ 
tive  functions  and  the  group  problem.  Recent  talks  have  been  given  on  this 
subject  [13].  Cutting  plane  generation  from  subadditive  functions  is  a  gen¬ 
eral  framework.  The  k- cuts  in  this  paper  can  be  viewed  as  a  particular  case 
of  this  framework;  the  subadditive  functions  for  £:-cuts  are  the  simplest  ones 
among  the  possible  subadditive  functions.  Nevertheless,  it  is  interesting  to 
study  this  particular  case,  because  it  comes  from  a  familiar  and  very  natural 
operation. 

In  a  section  of  their  book  [7]  Garfinkel  and  Nemhauser  discuss  the  pos¬ 
sibility  of  obtaining  stronger  cuts  than  the  Gomory  fractional  cuts  by  first 
multiplying  a  tableau  row  by  a  nonzero  factor  before  generating  the  fractional 
cut.  The  potential  multipliers  could  be  integral  or  non-integral  numbers.  For 
integral  multipliers,  the  larger  the  fractional  part  of  the  right-hand-side  of  the 
resulting  tableau  row,  the  stronger  the  resulting  cut  becomes  in  most  of  the 
cases.  For  non-integral  multipliers,  the  situation  turns  out  to  be  complicated. 

Dawande  and  Pulleyblank  [6]  considered  the  cone  formed  at  the  LP  opti¬ 
mum  by  the  normals  of  the  tight  constraints  at  the  optimum.  The  direction 
vector  (the  normal  to  the  cut  hyperplane)  of  the  Gomory  fractional  cut  lies 
in  this  cone.  The  direction  vector  is  integral.  In  general,  there  are  several 
integral  vectors  in  the  cone.  After  multiplying  the  tableau  rows  by  inte¬ 
gral  multipliers,  precisely  which  integral  vectors  in  the  cone  are  the  direction 
vectors  of  the  cuts?  Dawande  and  Pulleyblank  have  done  some  work  on 
characterizing  these  directions  in  order  to  better  understand  the  Gomory 
fractional  cuts  from  a  geometrical  point  of  view. 

Giinliik  and  Pochet  [17]  described  a  simple  cutting-plane  algorithm  that 
uses  a  “mixing  procedure”  to  generate  valid  inequalities.  First,  they  generate 
10  base  inequalities  from  each  tableau  row  by  multiplying  the  row  with  the 
multipliers  <5  =  1,2,  3,  4,  5  and  5  =  —  1,  —2,  —3,  —4,  —5.  And  then  they  apply 
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their  “strengthening  procedure”  to  the  collection  of  the  base  inequalities  to 
obtain  cuts  and  improve  the  relaxation  bound  for  general  integer  programs. 

In  this  paper  we  focus  on  a  variation  of  Gomory  mixed  integer  cuts  for 
pure  integer  programs.  Rather  than  generating  Gomory  mixed  integer  cuts 
from  the  LP  tableau  rows  directly,  we  first  multiply  the  tableau  rows  by  an 
integer  k  and  then  generate  the  cuts  from  the  new  rows  according  to  the 
Gomory  mixed  integer  procedure.  We  call  them  k-cuts  in  this  paper.  A 
comparison  of  £:-cuts  (k  >  2)  with  the  usual  Gomory  mixed  integer  cuts 
k  =  1)  is  given  in  Section  3.  It  is  proved  in  the  pure  case  that  with  exactly 
50%  probability  the  ft-cuts  perform  variable-wise  better  than  the  Gomory 
mixed  integer  cuts.  These  fc-cuts  provide  a  variety  of  cuts  different  from 
Gomory  mixed  integer  cuts  and  can  cut  off  the  relaxation  set  from  different 
angles.  Computational  results  on  knapsack  problems  are  given  in  Section  4 
to  illustrate  that  the  £:-cuts  have  the  same  power  as  Gomory  mixed  integer 
cuts  in  terms  of  improving  the  relaxation  bound.  If  we  impose  all  the  ft-cuts 
(1  <  k  <  M )  as  constraints,  we  see  empirically  that  they  are  effective  as 
a  whole  in  restricting  the  feasible  set  of  the  LP  relaxation  of  the  knapsack 
problem  when  M  becomes  larger.  Surprisingly,  the  same  does  not  appear  to 
be  true  for  integer  programs  with  many  constraints,  as  illustrated  at  the  end 
of  Section  4.  An  extension  of  ft-cuts  to  mixed  integer  programs  is  discussed 
in  the  last  section. 

2.  Gomory  mixed  integer  cuts  and  k- cuts 

In  this  section,  we  start  with  the  definition  and  properties  of  ft-cuts,  and  we 
give  two  examples.  Finally,  we  introduce  a  graphic  tool  that  is  helpful  when 
we  compare  the  cuts  in  the  next  section. 

Consider  the  pure  integer  program  (IP) 

Min  cx 
s.t.  Ax  =  b, 
x  G  Z", 

where  A  is  a  rational  matrix,  and  c  and  b  are  rational  vectors. 

After  solving  the  LP  relaxation  of  (IP),  we  have  the  following  basic  solu¬ 
tion 

Xi  =  ai o  +  Ejej  aij(-Xij)  for  all  i  G  /,  (1) 

where  I  and  J  denote  the  sets  of  basic  and  non-basic  variables. 

Write  aij  =  [aij J  +  fij  (i  G  I  and  j  G  J)  and  ai0  =  [ai0\  +  fi0  (i  G  /). 
Assume  /,0  >  0  for  some  i.  The  Gomory  fractional  cut  (GFC)  from  that  LP 
tableau  row  is 
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Ejej  fijXj  >  fio-  (2) 

The  Gomory  mixed  integer  cut  (GMIC)  from  the  LP  tableau  row  is: 

fijXj  +  l- fio  ~  fij)xj  —  fio ■  (3) 

Rather  than  directly  generating  the  GMIC  from  the  LP  tableau  row,  we 
could  first  multiply  both  sides  of  the  row  by  some  nonzero  integer  k  ^  1 
and  then  apply  the  Gomory  mixed  integer  procedure  as  above  to  obtain  the 
corresponding  inequality: 

^3xi  +  I ^ (j >/«)(-*-  -  fij)xi  -  /«,  (4) 

where  fi0  =  ai0  -  [ai0\  >  0,  fij  =  -  [a^J  >  0  (j  G  J),  ai0  =  kai0  and 

dij  =  kdij  ( j  e  J).  We  call  the  inequality  in  (4)  a  k-cut  generated  from  the 
LP  tableau  row. 

LEMMA  1.  The  k-cut  is  a  valid  inequality  for  (IP). 

PROOF.  Multiplying  both  sides  of  the  LP  tableau  row  by  a  nonzero 
integer  k  ^  0,  we  have  kx ,  =  a,0  +  J2jeJ  dij(~xij)-  Replacing  kx j  by  a  new 
integer  variable  x'i:  we  get  x\  =  a,0  +  J2jeJ^ij(~xij)i  which  is  an  equation 
similar  to  (1).  The  standard  derivation  of  GMIC  applied  to  this  equation 
shows  that  (4)  is  valid  for  IP.  □ 

LEMMA  2.  The  k-cut  is  equivalent  to  the  (-k)-cut. 

PROOF.  It  is  sufficient  to  show  that  the  GMIC  is  equivalent  to  the  (— 1)- 
cut.  The  GMIC  from  the  LP  tableau  row  is  shown  in  (3),  which  is 

Y  faxi  +  Y  fiOxj  +  73V  ^  “  fij)X3  -  /«■ 

After  multiplying  —1  to  the  tableau  row,  the  fractional  parts  of  the  coeffi¬ 
cients  in  the  new  row  are  just  1  —  fij  and  1  —  /,0.  Therefore,  the  (— l)-cut 
is 


Y  C1  -  fij)xj  +  Y  (4  “  fio)xj  +  1  f/t0  Y  fijxj  >  1  -  fio- 

We  now  see  that  there  exists  just  a  positive  factor  between  the  GMIC 
and  the  (— l)-cut.  □ 

Now  let  us  try  to  get  some  sense  of  the  variety  and  efficiency  of  the  fc-cuts 
by  analyzing  two  small  examples. 
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EXAMPLE  1  (general  integer  case) 

Consider  the  following  integer  programming  problem: 

Max  ICte  +  13 y 
s.t.  lOx  +  14 y  <  43, 
x,  y  e  Z+. 

It  is  equivalent  to  the  following  IP  after  introducing  an  integer  slack  z: 

Min  y  +  z 

s.t.  x  +  lAy  +  O.lz  =  4.3, 
x,  y,  z  e  Z+. 

Obviously,  the  equality  x  +  1.4y  +  0.1  z  =  4.3  is  a  row  of  the  corresponding 
optimal  LP  tableau.  The  LP  optimum  is  ( x ,  y,  z )  =  (4.3,  0,  0). 

All  the  possible  k- cuts  generated  from  this  tableau  row  are  summarized 
in  the  following  table  and  Figure  1. 


GMIC  and  9-cut 

f  y  +  2z  >  6 

2-cut  and  8-cut 

3y  +  2z  >  6 

3-cut  and  7-cut 

2y  +  2>z  >  9 

4-cut  and  6-cut 

2y  +  3z  >  4 

5-cut 

z  >  1 

A  picture  of  these  £:-cuts  in  the  (y,  z)-plane  clearly  shows  that  all  the 
ft-cuts  (2  <  k  <  8)  either  dominate  the  GMIC  or  are  incomparable  with  the 
GMIC  in  this  example.  These  ft-cuts  effectively  shrink  the  feasible  region  of 
the  LP  relaxation  by  cutting  off  different  parts.  Worth  mentioning  is  that 
the  3-cut  is  a  facet  to  the  convex  hull  of  the  feasible  integer  points  in  the 
(y,  z)-plane.  Together  with  the  valid  inequalities  14 y  +  z  <  43  and  y  >  0, 
it  gives  a  complete  facet  description  of  the  convex  hull  of  the  feasible  (y,  z) 
integer  points. 

If  we  consider  the  following  IP  instead  of  the  above  one 

Min  y  +  z 

s.t.  x  +  OAy  +  O.lz  =  5.3, 
x,  y,  z  e  Z+, 

then  all  the  £:-cuts  generated  from  the  new  optimum  LP  tableau  row  are  the 
same  as  above.  However,  in  this  case  the  3-cut  and  the  5-cut,  together  with 
4  y  +  z  <  53  and  y  >  0,  give  a  complete  facet  description  of  the  corresponding 
convex  hull  of  the  feasible  (y,  z )  integer  points.  □ 
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Figure  1:  £:-cuts  in  (y,  z)-plane  for  Example  1 


EXAMPLE  2  (0-1  case) 

Consider  the  following  0-1  programming  problem: 

Max  xi  +  x2 
s.t.  10xi  +  2x2  <  11, 

2xi  +  10x2  <  11, 

Xi  <  1,  x2  <  1, 

xi,  x2  e  z+. 

After  introducing  two  integer  slacks  w  and  z,  we  get  the  equivalent  IP: 

Max  xi  +  x2 
s.t.  10xi  +  2x2  +  w  =  11, 

2a:  i  +  10x2  +  z  =  11, 

Xi  <  1,  x2  <  1, 

Xi,  x2  w,  zG  Z+. 

The  GMICs  generated  from  the  optimum  LP  tableau  rows  Xi+^w—^z  = 
and  x2  -  ^ w  +  are 

f  72Xl  +  120^2  <  132, 

1  120xi  +  72x2  <  132. 
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Figure  2:  9-cuts  vs  GMIC  in  (xi,  a;2)-plane  for  Example  2 


Their  respective  boundary  lines  pass  through  the  points  “(1,  |)  and  (|,  1)” 
and  “(1,  |)  and  (|,  1)”  and  they  intersect  at  the  point  ff )  in  the  (xi,x2)- 

space. 

The  9-cuts  generated  from  those  two  LP  tableau  rows  are 


X\  +  2x2  <  2, 

2xi  T  X2  T  2. 


The  corresponding  boundary  lines  pass  through  the  points  “(1,  |)  and  (0, 1)” 
and  “(1,0)  and  (|,  1)”  respectively  and  intersect  at  (§,§)• 

We  see  that  these  two  9-cuts  respectively  dominate  the  corresponding 
GMICs.  In  fact,  together  with  the  non-negativity  constraints  on  x\  and  x2, 
these  two  9-cuts  give  a  complete  description  of  the  elementary  closure  of  the 
LP  relaxation  under  the  lift-and-project  cuts.  □ 

3.  Comparing  GMICs  with  £>cuts 

In  this  section  we  are  going  to  discuss  the  relation  between  Gomory  mixed 
integer  cuts  and  fc-cuts  in  pure  integer  programs.  The  relation  builds  upon 
the  following  definition. 

DEFINITION  1.  Let  ax  >  (3  and  ax  >  j3  be  two  cuts  for  (IP),  where 
a,  a  €  Rf,  (3  >  0  and  j3  >  0.  The  cut  ax  >  (3  is  equally  strong  on  Xj  as  the 
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cut  ax  >  (3,  if  aj/ /3  =  aj//3.  The  cut  ax  >  (3  is  stronger  on  Xj  than  the  cut 
ax  >  (3,  if  aj/ (3  <  aj/(3.  The  cut  ax  >  j3  is  strictly  stronger  on  Xj  than  the 
cut  ax  >  (3,  if  aj/ (3  <  a.j/f3.  □ 

Obviously,  if  the  cut  ax  >  (3  is  stronger  on  Xj  than  the  cut  ax  >  f3  for 
all  1  <  j  <  n,  then  the  cut  ax  >  (3  is  stronger  than  the  cut  ax  >  f3,  or  in 
other  words,  ax  >  f3  dominates  ax  >  (3. 

Let  x0  =  a0  —  YjjejhjXj  be  an  LP  tableau  row,  where  aj  =  |_ aj\  +  fj, 
ho  =  L«0J  +  /o  and  0  <  /0  <  1.  Let  Y,jeJajxj  ^  P  represent  any  of  the 
three  cuts  generated  from  this  tableau  row:  GFC  as  described  in  (2),  GMIC 
as  described  in  (3)  and  the  £:-cut  as  described  in  (4).  Then  aj  >  0  (je  J) 
and  j3  >  0.  Given  0  <  /0  <  1,  aj/(3  is  in  each  case  only  a  function  of 
fj.  Actually,  the  function  aj/(3  for  GFC  is  a  linear  function  starting  from 
origin,  the  function  for  GMIC  can  be  plotted  as  a  two-slope  figure,  and  the 
function  for  a  k- cut  ( k  >  2)  is  a  piecewise  linear  function  with  a  regular  zig¬ 
zag  feature.  For  instance,  when  /0  =  0.15,  the  corresponding  functions  for 
GFC  and  GMIC  are  plotted  in  Figure  3.  The  function  ct,//3  corresponding  to 
the  5-cut  when  /0  =  0.15  is  shown  in  Figure  4.  This  kind  of  picture  allows  to 
compare  two  inequalities  by  plotting  the  two  functions  and  checking  where 
one  curve  is  below  the  other. 


Figure  3:  aj/ (3  as  a  function  of  fj  for  GFC  and  GMIC 
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Figure  5  illustrates  the  relation  among  the  functions  cxj/ f3  for  GFC,  GMIC 
and  a  5-cut  when  /0  =  0.15.  Since  k  =  5  and  /0  =  0.15,  we  have  k  < 

We  can  observe  from  the  figure:  (1)  The  GFC,  GMIC  and  the  5-cut  are 
equally  strong  on  Xj  if  0  <  fj  <  /0;  (2)  The  GMIC  and  the  5-cut  are  strictly 
stronger  on  Xj  than  GFC  if  /0  <  fj  <  1;  (3)  The  GMIC  is  strictly  stronger 
on  Xj  than  the  5-cut  in  some  of  the  sub-intervals  of  (/0, 1),  while  in  the  rest, 
except  a  few  isolated  points,  5-cut  is  strictly  stronger  on  Xj  than  GMIC. 

The  observation  from  the  example  in  Figure  5  can  be  extended  to  any 
0  <  /0  <  0.5  and  k  <  |_^_|-  The  comparison  is  summarized  in  Table  1,  where 
we  use  the  numbers  1,  2  and  3  to  rank  the  strength  of  GFC,  GMIC  and  ft-cut 
on  Xj  in  a  decreasing  order. 

It  is  easy  to  sum  up  from  Table  1  the  total  length  of  the  sub-intervals 
where  the  £:-cut  is  strictly  stronger  on  Xi  than  the  GMIC.  This  length  is  equal 
to  This  exactly  equals  the  total  length  of  the  sub-intervals  where  the 

GMIC  is  strictly  stronger  on  Xj  than  the  £:-cut. 


fj  value 

A:-cut 

GMIC 

GFC 

[0,/o] 

1 

1 

1 

(/o,  i: 

1 

2 

3 

r  i  j+i'i 

U’  k  /’ 

1  <  i  <  k  -  1 

[f,  f  +  /o] 

%  _i_  k— i  f  \ 

1*’  k  "T  k  JO) 

1 

2 

t  i  k—i  £ 

k  k  ^0 

1 

1 

(1  +  nrfoi  i  +  fo] 

2 

1 

(i  +  /o, 

(fc  +  fo,  fc_!  1  k—1  fo) 

2 

1 

1  |  k—i— 1 x  ^ 

k- 1  '  k- 1  JO 

1 

1 

7  *  i  k—i— l £  j+i \ 
Vfe-1  “r  k- 1  JOi  k  > 

1 

2 

Table  1:  comparison  among  cuts  when  0  <  /0  <  0.5  and  1  <  k  <  |_^J 

THEOREM  1.  Let  x0  =  o0  —  Jf.jejdjXj  be  a  LP  tableau  row,  where 
aj  =  [ hj\  +  fj,  o0  =  L«oJ  +  fo  and  0  <  /o  <  0.5.  Then  we  have 

(i)  When  1  <  k  <  ^ne  k-cut  generated  from  the  row  is  stronger  than 
the  GFC  generated  from  the  row. 

(ii)  Assume  that  f  is  uniformly  distributed  in  [0, 1).  When  2  <  k  < 

the  probability  that  the  k-cut  generated  from  the  row  is  strictly  stronger  on 
Xi  than  the  GMIC  generated  from  the  row  is  which  is  equal  to  the 

probability  that  the  GMIC  generated  from  the  row  is  strictly  stronger  on  Xj 
than  the  k-cut  generated  from  the  row.  □ 
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When  we  compare  the  £:-cuts  (2  <k<  L^J)  as  a  group  with  the  GMIC, 
we  have  the  following  result. 

THEOREM  2.  A  ssume  that  <  fa  <  for  some  integer  K  >  2 
and  that  fa  is  uniformly  distributed  in  [0,1).  The  probability  that  some  k- 
cut  generated  from  the  row  (2  <  k  <  K)  is  strictly  stronger  on  Xj  than  the 
GMIC  generated  from  the  row  is  while  the  probability  that  the 

GMIC  generated  from  the  row  is  strictly  stronger  on  Xj  than  all  the  k-cuts 
generated  from  the  row  (2  <  k  <  K)  is  — j 

PROOF.  From  Table  1  we  see  that  the  GMIC  and  the  fc-cut  (2  <  k  <  K ) 
are  equally  strong  on  Xj  when  0  <  fj  <  fa- 

Table  1  shows  that  the  2-cut  is  strictly  stronger  on  Xj  than  the  GMIC 
when  fa  <  fj  <  In  the  case  of  k  =  2  and  i  =  1,  the  2-cut  is  strictly 
stronger  on  Xj  than  the  GMIC  when  \  <  fj  < 

Let  3  <  p  <  K.  In  the  case  of  k  =  p  and  i  =  p  —  2,  the  p-cut  is  strictly 
stronger  on  Xj  than  the  GMIC  when  p~pff/°  <  fj  <  while  in  the  case  of 
k  =  p  and  i  —  p  —  1  the  p-cut  is  strictly  stronger  on  Xj  than  the  GMIC  when 
<  fj  <  p~ lpf° .  So  the  p-cut  is  strictly  stronger  on  Xj  than  the  GMIC 
when 

p— 1  J J  p 

If  we  put  together  all  the  intervals  (/0,  |),  [|,  and  ,  p  l^° ) 

(3  <  p  <  K),  we  see  that,  for  any  given  value  of  fj  €  (/0,  )  except  on 

K  —  2  isolated  points  p~1^0  (2  <  p  <  K  —  1),  there  must  exist  some  ft-cut 
(2  <  k  <  K)  such  that  the  k- cut  is  strictly  stronger  on  Xj  than  the  GMIC. 
The  length  of  the  interval  (/0,  )  is  . 

It  is  not  difficult  to  see  in  the  case  of  k  =  p  and  i=p— l(2<p<K) 
that  the  GMIC  is  strictly  stronger  on  Xj  than  the  p-cut  when  fj  e  (p~p+^°,  1)- 
Therefore,  the  GMIC  is  strictly  stronger  on  Xj  than  all  the  A;-cuts  (2  <  k  < 
K )  when  fj  €  (ffizThio  ^  l).  The  length  of  the  interval  ( ,  1)  is  □ 

Now  let  us  compare  the  GMIC  with  a  £:-cut  in  the  case  kf0  >  1. 

In  Figure  6  (/0  =  0.375)  and  Figure  7  (/0  =  0.875),  we  can  respectively 
observe  the  relation  between  the  functions  otj/fi  of  the  GMIC  and  the  4- 
cut.  We  have  \kf0]  =  [4  x  0.375]  =  [1.5]  =  2  <  4  =  k  in  Figure  6  and 
pfe/ol  =  [4  x  0.875]  =  [3.5]  =  4  =  k  in  Figure  7.  In  fact,  these  two  figures 
effectively  represent  all  the  possible  combinations  of  k  and  fa  when  kfa  >  1. 
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Table  2:  comparison  between  ft-cut  and  GMIC  when 
0  <  /  <  kf0  <  l  +  1  <  k 


Table  3:  comparison  between  k-cut  and  GMIC  when 
0  <  l  <  kf0  <  l  +  1  =  k 


The  observations  in  Figure  6  and  Figure  7  can  be  quantitatively  summa¬ 
rized  in  Table  2  and  Table  3  respectively  in  the  general  cases  0  <  l  <  kf0  < 
l  +  1  <  k  and  0  <  l  <  kf0  <  l  +  1  =  k  for  some  integer  l.  Again  we  use 
the  numbers  1  and  2  in  the  tables  to  rank  the  strength  of  the  cuts  on  Xj 
in  a  decreasing  order.  Note  that  the  case  listed  as  fj  =  1  in  Table  2  and 
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Table  3  corresponds  to  the  actual  case  fj  =  0.  Now  Theorem  3  follows  from 
observing  the  figures  and  tables. 

THEOREM  3.  Assume  that  fa  is  uniformly  distributed  in  [0,1).  Let 
k  >  |~y-~|  and  kf0  be  non-integral. 

(i) Ifk  \k ,/o  |  j  then  the  probability  that  the  k  cut  generated  from  the 
tableau  row  is  strictly  stronger  on  Xi  than  the  GMIC  generated  from  the  row 
is  which  is  equal  to  the  probability  that  the  GMIC  generated  from  the  row 
is  strictly  stronger  on  x j  than  the  k-cut  generated  from  the  row. 

(ii)  If  k  =  ffe/ol ,  then  the  probability  that  the  k-cut  generated  from  the 

tableau  row  is  strictly  stronger  on  Xi  than  the  GMIC  generated  from  the  row 
is  which  is  equal  to  the  probability  that  the  GMIC  generated  from  the  row 
is  strictly  stronger  on  Xj  than  the  k-cut  generated  from  the  row.  □ 

In  the  following,  we  compare  the  GMIC  to  a  group  of  &:-cuts  when  /0  > 

THEOREM  4.  Assume  fo  >  1  and  q  <  mfo  <  q  +  1  for  some  integer 
q  >  1  and  integer  m  >  2.  Then 

(i)  When  fj  €  (0,  the  GMIC  is  strictly  stronger  on  Xj  than  any  k-cut 
(2  <  k  <  m). 

(ii)  When  fj  £  fo],  there  exists  some  k-cut  (2  <  k  <  m)  such 

that  the  k-cut  is  stronger  on  Xj  than  the  GMIC,  and  the  “strictly  stronger” 
property  holds  for  all  fj  £  (^p,  fo]  except  a  finite  number  of  isolated  points. 

(Hi)  If  f0  >  then  the  2-cut  and  the  GMIC  are  equally  strong  on  Xj 
when  fj  £  [fo,  1).  If  fo  =  then  the  3-cut  is  strictly  stronger  on  Xj  than  the 
GMIC  when  fj  £  (|,  f). 

PROOF,  (i)  In  the  case  of  *  =  0  in  Table  2  and  Table  3,  the  GMIC  is 
strictly  stronger  on  Xj  than  any  k- cut  (2  <  k  <  m)  when  fj  £  (0,  ^[p]. 

(ii)  Let  k'  be  the  largest  integer  between  2  and  m  such  that  \k' fo]  =  k' . 
First,  let  us  look  at  Table  2  for  the  comparison,  where  we  only  consider  the 
integer  value  of  k  in  [k'  +  1,  m\.  We  decrease  the  value  of  k  from  m  to  k'  +  1 
in  a  way  to  ensure  that  the  integer  l  decreases  one  at  a  time.  Each  time 
when  l  has  a  lower  value,  we  choose  i  =  0  and  1  to  identify  two  sub-intervals, 
in  which  the  ft-cut  is  strictly  stronger  on  Xj  than  the  GMIC.  If  we  connect 
together  all  these  sub-intervals,  then  we  get  an  interval  (^p,  pAp].  From  how 
we  construct  this  interval,  we  know  that,  when  fj  £  (^p,  p/zp]>  there  exists 
some  k- cut  (k'  +  1  <  k  <  m)  such  that  the  &:-cut  is  stronger  on  Xj  than 
the  GMIC,  and  the  “strictly  stronger”  property  holds  for  all  fj  £  (^p,  pAp] 
except  a  few  isolated  points. 

Now  let  us  compare  the  fc-cuts  with  GMIC  in  Table  3,  assuming  that  the 
integer  value  of  k  is  taken  in  [2,  k'].  Similarly  as  how  we  identify  and  connect 


14 


the  sub-intervals  above,  we  decrease  the  value  of  k  once  a  time,  meanwhile 
the  value  of  l  decreases  once  a  time  as  well.  Each  time  we  choose  i  =  0  and 
1  to  figure  out  two  sub-intervals,  in  which  the  k-cut  is  strictly  stronger  on 
Xj  than  the  GMIC.  Connecting  together  all  these  sub-intervals,  we  obtain 
an  interval  (jj,  |].  Also  we  see  in  Table  3  that  2-cut  is  strictly  stronger  on 
Xj  than  the  GMIC  when  fj  e  (|,/0).  Therefore,  when  fj  e  (|j,/o],  there 
exists  some  &-cut  (2  <  k  <  k')  such  that  the  &-cut  is  stronger  on  Xj  than  the 
GMIC,  and  the  “strictly  stronger”  property  holds  for  all  fj  €  except 

a  few  isolated  points. 

(iii)  In  the  case  of  /0  >  |,  Table  3  tells  us  that  the  2-cut  and  the  GMIC 
are  equally  strong  on  Xj  when  fj  e  [/0, 1). 

In  the  case  of  /0  =  |,  the  3-cut  is  strictly  stronger  on  Xj  than  the  GMIC 
when  fj  G  (|,  |),  which  can  be  seen  if  we  look  at  the  case  of  k  =  3,  l  =  1 
and  i  =  2  in  Table  2.  □ 

4.  Computational  results 

In  this  section,  we  test  computationally  how  the  &-cuts  compare  to  GMIC  for 
0-1  and  bounded  knapsack  problems  and  for  integer  programs  with  multiple 
constraints. 

The  knapsack  problems  have  the  following  form: 

Max  YJj=iPjXj 
S.t.  Yfj= 1  wjXj  <  C, 

0  <  Xj  <  bj  and  integer,  1  <  j  <  n, 

where  pj,  Wj  and  c  are  positive  numbers  and  bj  is  a  positive  integer.  The  0-1 
knapsack  problems  have  bj  =  1  for  all  1  <  j  <  n. 

For  the  0-1  knapsack  problems,  we  generate  both  pj  and  Wj  randomly 
uniformly  in  [1,1000]  and  c  =  0.5X)"=i  Wj,  following  Martello  and  Toth  [19]. 
When  we  generate  the  bounded  knapsack  problems,  pj  and  Wj  are  set  to  be 
uniformly  random  in  [1,1000],  the  values  bj  are  uniformly  random  in  [5,10], 
and  c  is  set  to  0.5  Yfj=\  bjWj.  The  number  n  of  variables  ranges  from  10  to 
10000,  which  represents  the  size  of  the  problem. 

In  Table  4  and  Table  5,  we  compare  the  gap  closed  after  adding  one  £:-cut 
to  the  LP  relaxation  of  the  knapsack  problems,  where  k  ranges  from  1  to 

5.  The  percentage  figures  in  the  tables  represent  the  gap  between  the  LP 
optimum  and  the  IP  optimum  closed  by  adding  the  ft-cuts.  The  asterisks 
V  in  the  upper  right  corner  of  some  percentage  values  indicate  those  values 
that  are  the  largest  among  all  the  percentage  values  in  their  row.  These  two 
tables  show  that,  for  both  the  0-1  and  the  bounded  cases,  each  £:-cut  has 
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the  same  chance  of  being  the  best  cut,  over  all  values  of  k.  The  GMICs 
do  not  show  any  extraordinary  advantage  over  other  ft-cuts.  This  fact  is  in 
accordance  with  the  results  in  Theorem  1  and  Theorem  3  in  Section  3. 


Problem 

set 

Number  of 
variables 

gap  closed 
by  1-cut 

gap  closed 
by  2-cut 

gap  closed 
by  3- cut 

gap  closed 
by  4-cut 

gap  closed 
by  5-cut 

problem  1 

10 

76.31% 

56.79% 

92.22% 

32.39% 

problem  2 

10 

|  56.88% 

29.82% 

17.64% 

6.18% 

problem  3 

100 

J  5.94% 

8.32% 

3.50% 

3.33% 

problem  4 

100 

14.87% 

15.66% 

28.72% 

23.21% 

problem  5 

1000 

12.25% 

6.06% 

5.20% 

25.74% 

problem  6 

1000 

13.15% 

memmm 

9.51% 

0.22% 

5.20% 

problem  7 

5000 

40.74% 

34.26% 

46.30% 

55.56% 

problem  8 

5000 

9.01% 

9.91% 

8.11% 

11.71% 

problem  9 

10000 

2.94% 

4.41% 

5.88% 

problem  10 

10000 

7.22% 

21.65% 

;  8.25% 

7.22% 

Table  4:  individual  ft-cuts  in  0-1  knapsack  problems 


Problem 

set 

Number  of 
variables 

gap  closed 
by  1-cut 

gap  closed 
by  2-cut 

gap  closed 
by  3- cut 

gap  closed 
by  4-cut 

gap  closed 
by  5-cut 

problem  1 

10 

84.30% 

70.17% 

60.92% 

52.18% 

problem  2 

10 

19.73% 

21.51% 

32.19% 

6.48% 

problem  3 

100 

13.94% 

13.94% 

21.91% 

29.51% 

problem  4 

100 

49.90% 

27.29% 

42.70% 

!  43.18%  j 

problem  5 

1000 

18.01% 

problem  6 

1000 

4.11% 

2.57% 

i  1.80% 

3.34% 

problem  7 

5000 

8.26% 

18.18% 

14.05% 

14.88% 

problem  8 

5000 

1.42% 

0.71% 

0.71% 

0.71% 

problem  9 

10000 

9.19% 

5.95% 

5.41% 

6.49% 

■IfcbtaAl 

problem  10 

10000 

2.68% 

4.03% 

8.05% 

;  2.01%  i 

Table  5:  individual  ft-cuts  in  bounded  knapsack  problems 


Problem 

set 

Number  of 
variables 

gap  closed 
by  1-cut 

gap  closed 
by  [l,10]-cuts 

gap  closed 
by  [l,50]-cuts 

problem  1 

10 

53.98% 

92.57% 

100% 

problem  2 

10 

94.05% 

100% 

100% 

problem  3 

50 

19.19% 

37.58% 

40.24% 

problem  4 

50 

8.11  % 

35.01% 

35.64% 

problem  5 

100 

30.78% 

77.04% 

86.18% 

problem  6 

100 

53.09% 

100% 

100% 

problem  7 

500 

18.39% 

25.06% 

52.18% 

problem  8 

500 

16.39% 

49.76% 

55.06% 

problem  9 

1000 

6.79% 

43.82% 

54.58% 

problem  10 

1000 

92.16% 

100% 

100% 

problem  11 

5000 

15.77% 

45.00% 

57.31% 

problem  12 

5000 

15.27% 

39.69% 

74.05% 

problem  13 

10000 

3.70% 

17.78% 

24.44% 

problem  14 

10000 

6.35% 

14.81% 

23.28% 

Table  6:  accumulated  £-cuts  in  0-1  knapsack  problems 
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Problem 

set 

Number  of 
variables 

gap  closed 
by  1-cut 

gap  closed 
by  [l,10]-cuts 

gap  closed 
by  [l,50]-cuts 

problem  1 

10 

66.95% 

100% 

100% 

problem  2 

10 

74.29% 

100% 

100% 

problem  3 

50 

41.87% 

64.95% 

77.30% 

problem  4 

50 

6.45% 

66.51% 

70.99% 

problem  5 

100 

39.62% 

66.00% 

69.04% 

problem  6 

100 

20.38% 

58.16% 

94.71% 

problem  7 

500 

15.25% 

43.17% 

71.29% 

problem  8 

500 

2.22% 

2.22% 

66.37% 

problem  9 

1000 

6.30% 

20.82% 

39.45% 

problem  10 

1000 

6.46% 

39.38% 

81.88% 

problem  11 

5000 

8.62% 

17.24% 

25.00% 

problem  12 

5000 

11.56% 

37.41% 

41.50% 

problem  13 

10000 

2.80% 

6.54% 

17.76% 

problem  14 

10000 

1.67% 

16.67% 

36.67% 

Table  7:  accumulated  &-cuts  in  bounded  knapsack  problems 


Problems  with 

Number  of 

gap  closed 

gap  closed 

gap  closed 

500  0-1  var. 

constraints 

by  1-cut 

by  [1, 10]-cuts 

by  [l,50]-cuts 

problem  1 

5  constr. 

12.53% 

15.56% 

15.62% 

problem  2 

10  constr. 

1.13% 

2.27% 

2.57% 

problem  3 

50  constr. 

4.04% 

4.16% 

4.16% 

problem  4 

100  constr. 

2.32% 

2.32% 

2.32% 

probelms  with 

Number  of 

gap  closed 

gap  closed 

gap  closed 

500  bd  var. 

constraints 

by  1-cut 

by  [1, 10]-cuts 

by  [l,50]-cuts 

problem  5 

5  constr. 

2.42% 

5.33% 

5.46% 

problem  6 

10  constr. 

5.09% 

5.09% 

5.09% 

problem  7 

50  constr. 

0.75% 

0.79% 

0.79% 

problem  8 

100  constr. 

1.88% 

1.92% 

1.92% 

Table  8:  accumulated  fc-cuts  in  integer  programming  problems 

In  Table  6  and  Table  7,  we  compare  the  cuts  cumulatively:  the  third 
column  corresponds  to  the  case  of  just  adding  the  GMIC,  the  next  column 
corresponds  to  adding  all  the  &-cuts  with  1  <  k  <  10,  and  the  last  column 
corresponds  to  adding  all  the  £:-cuts  with  1  <  k  <  50.  The  data  imply  that 
adding  several  £:-cuts  can  close  the  gap  effectively  and  works  better  than  any 
single  fc-cut.  In  some  cases,  the  improvement  is  very  significant.  The  results 
in  Theorem  2  and  Theorem  4  are  partly  reflected  in  these  tables. 

In  Table  8  we  consider  applying  k- cuts  to  the  0-1  and  bounded  integer 
programming  problems  with  more  than  one  constraint.  Each  problem  has 
500  variables  and  from  5  to  100  constraints.  The  objective  function  and  each 
constraint  are  randomly  generated  as  in  the  problems  in  Tables  4-7.  We 
generate  &:-cuts  from  all  the  rows  of  the  tableau  in  which  the  basic  variable 
has  a  fractional  value.  The  results  of  Table  8  show  that  the  effect  of  the 
&:-cuts  tails  off  dramatically  as  the  number  of  constraints  increases.  This 


17 


is  rather  disappointing.  Even  more  disappointing  is  the  observation  that 
adding  groups  of  k- cuts  provides  only  marginal  improvements  over  adding 
the  GMICs  alone.  We  do  not  have  any  good  theoretical  explanation  for  this. 

Obviously,  the  problems  with  more  than  one  constraint  are  more  mean¬ 
ingful  from  a  practical  point  of  view.  However,  in  this  case,  adding  the  £:-cuts 
for  all  k  <  50  is  not  efficient  as  clearly  illustrated  in  Table  8.  A  further  in¬ 
vestigation  of  &-cuts  might  still  be  interesting,  especially  for  problems  with 
general  integer  variables,  since  there  are  not  many  powerful  cuts  available 
for  the  general  integer  case.  To  investigate  the  usefulness  of  using  fc-cuts  in 
practice,  one  should  solve  problems  to  optimality,  say  in  combination  with 
the  branch-and-bound  method.  Our  purpose  in  this  section  was  much  more 
limited.  We  only  wanted  to  illustrate  the  theorems  of  Section  3. 

5.  Extension  to  the  mixed  case 

The  mixed  integer  programming  (MIP)  problems  have  the  standard  form 

Min  cx  +  dy 
s.t.  Ax  +  By  =  b, 
xeZ^,  ye  Rp+, 

where  A  and  B  are  rational  matrices,  and  c,  d  and  b  are  rational  vectors. 
Solving  the  LP  relaxation  of  the  MIP  gives  tableau  rows  of  the  form 

x0  =  a0  +  J2jeJ  ®j(—xj)  +  J2jes  9j(—yj),  (5) 

where  /0  =  a0  —  |_GoJ  >  0,  J  and  S  denote  the  sets  of  the  non-basic  integer 
variables  and  the  non-basic  continuous  variables. 

Using  the  notation  from  Section  2,  the  GMIC  is,  in  the  mixed  case: 

E  fjx3  +  UzV  E  i1  -  fj)xj  +  E  9jVj  -  zzz  E  9jVj  >  fo-  (6) 

r-fj<fo  1  10  j:gj> 0  1  j:9i< 0 


Note  that  the  difference  between  the  GMICs  in  (3)  and  (6)  is  that  the  cut  in 
(6)  contains  terms  for  the  continuous  variables. 

Applying  the  same  proof  procedure  as  that  for  Lemma  1  shows  that  the 
following  cutting  plane  is  valid  for  MIP: 

E  fix3  +  TzV  E  (!  -  fj)xJ  +  E  k9jVj  -  tzV  E  k9j9j  >  fo,  (7) 

j:fj<fo  1  j-Jj>f0  j:gj>0  1  JO  j:gj< 0 

which  is  called  the  k-cut  in  the  mixed  case,  or  the  ft-cut  in  brief. 
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Comparing  the  cuts  in  (6)  and  (7),  we  can  see  that,  in  the  mixed  case,  the 
fc-cut  (k  >  2)  is  weaker  in  the  continuous  variables  than  the  GMIC.  Similar  to 
the  picture  introduced  in  Section  3  for  an  integer  variable,  one  can  construct 
a  picture  for  a  continuous  variable  yj.  Let  /3  =  /0  and  aj  =  kgj  when  gj  >  0, 
jzj^kgj  when  gj  <  0.  Then  the  function  atj/fi  is  a  linear  function  of  \gj\  with 
nonnegative  slope  and  the  slope  increases  as  k  increases,  which  weakens  the 
cut  on  yj.  Hence  it  seems  logical  that  increasing  k  too  much  leads  to  bad 
inequalities. 

We  observed  in  practice  that  fc-cuts  (k  >  2)  in  the  mixed  case  are  often 
not  as  good  as  the  GMICs.  So  there  exists  an  intrinsic  difficulty  in  extending 
the  results  of  the  pure  IP  case  to  the  MIP  case.  One  might  still  be  able  to 
use  them  conditionally.  For  example,  in  the  case  where  /0  <  \  and  gj  >  0  for 
all  j  (or  /0  ^  and  gj  ^  0  for  all  j  'J ,  the  results  of  Theorem  1  and  Theorem 

2  can  be  applied  to  obtain  k- cuts  (2  <  k  <  |_^_|)  (or  2  <  k  <  Lpr^J)  that 
are  as  powerful  as  the  GMIC  and  provide  a  variety  of  cutting  planes. 

The  following  two  tables  illustrate  the  deterioration  of  £:-cuts  in  the  mixed 
case  as  k  increases.  In  the  knapsack  problems,  the  slack  variables  are  consid¬ 
ered  to  be  continuous.  The  corresponding  GMIC  and  &-cuts  are  generated 
in  the  forms  of  (6)  and  (7).  The  quality  of  the  cuts  is  measured  by  the  gap 
closed  after  adding  one  cut.  We  can  see  from  Table  9  and  Table  10  that 
the  gap  closed  tends  to  decrease  as  k  increases,  which  implies  that  the  £:-cut 
deteriorates  as  k  becomes  larger. 


problem 

set 

gap  closed 
by  1-cut 

gap  closed 
by  10-cut 

gap  closed 
by  100-cut 

gap  closed 
by  1000-cut 

gap  closed 
by  10000-cut 

problem  1 

10.50% 

13.22% 

3.17% 

0.87% 

0.06% 

problem  2 

5.94% 

2.47% 

1.89% 

0.77% 

0.03% 

problem  3 

5.60% 

0.27% 

7.96% 

0.27% 

0.02% 

problem  4 

30.31% 

24.19% 

5.83% 

0.09% 

0.09% 

problem  5 

51.93% 

16.76% 

6.86% 

0.55% 

0.08% 

Table  9:  individual  fc-cuts  in  0-1  knapsack  problems 


problem 

set 

gap  closed 
by  1-cut 

gap  closed 
by  10-cut 

gap  closed 
by  100-cut 

gap  closed 
by  1000-cut 

gap  closed 
by  10000-cut 

problem  1 

39.62% 

25.72% 

5.32% 

0.43% 

0.03% 

problem  2 

19.69% 

44.52% 

11.26% 

3.52% 

0.38% 

problem  3 

11.78% 

44.48% 

7.58% 

0.40% 

0.00% 

problem  4 

13.02% 

0.32% 

2.34% 

0.98% 

0.06% 

problem  5 

42.30% 

8.46% 

18.93% 

2.11% 

0.15% 

Table  10:  individual  £:-cuts  in  bounded  knapsack  problems 
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A  further  extension  of  GMIC  combines  several  rows  from  the  LP  tableau 
with  integer  multipliers.  This  idea  has  attracted  attention  recently  [1,  4,  17, 
24],  This  topic  is  currently  under  further  investigation. 
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